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Extensions of average Hamiltonian theory to quantum com- 
putation permit the design of arbitrary Hamiltonians, allow- 
ing rotations throughout a large Hilbert space. In this way, 
the kinematics and dynamics of any quantum system may be 
simulated by a quantum computer. A basis mapping between 
the systems dictates the average Hamiltonian in the quan- 
tum computer needed to implement the desired Hamiltonian 
in the simulated system. The flexibility of the procedure is 
illustrated with NMR on 13 C labelled Alanine by creating the 
non-physical Hamiltonian a z a z a z corresponding to a three 
body interaction. 

PACS numbers: 03.67.-a,76.60.-k 



I. INTRODUCTION 

In the early 1980's, researchers such as Benioff, Ben- 
nett, Deutsch, Feynman and Landauer |]-|5| studied the 
possibility of performing computations using the princi- 
ples of quantum mechanics, and conjectured that a ma- 
chine based on these principles might be able to solve 
certain types of problems more efficiently than can be 
done on a conventional von Neumann computer. Later, 
Lloyd proposed that such a quantum computer might be 
built from an array of coupled two-state quantum sys- 
tems, each of which can store one quantum bit, or qubit, 
of information Shortly thereafter, Shor proved that 
a quantum computer would be capable of factorizing in- 
tegers in polynomial time, thereby showing that the ex- 
ponential number of degrees of freedom accessible to a 
quantum computer does indeed enable it to solve some 
problems more efficiently than is believed possible on a 
conventional machine H], 

An essential feature of a universal quantum computer 
is the ability to transform efficiently any initial state 
vector to any other state vector within a large Hilbert 
space. Such operations can be thought of as using a 
real-number (continuous) rather than a Turing-machine 
(digital) model of computation |)|. Algorithms can be 
tailored to take advantage of this as well as the paral- 
lelism from quantum superposition. Since any quantum 
system governed by an arbitrary Hamiltonian may be 
described by the paths taken by a set of basis vectors, a 



quantum information processor should be able to simu- 
late the evolution of any smaller quantum system up to 
a specified time point. 

Accordingly, one of the first proposed applications of 
quantum computers was quantum simulation: using a 
quantum mechanical computer to simulate another quan- 
tum mechanical system. Feynman's original proposal in 
1982 H that there might be a universal and efficient 
quantum simulator of physical systems was recently val- 
idated, in general terms, by Lloyd [^0[ . Algorithms for 
specific classes of quantum systems have also been pro- 
posed 1 11 - 14 1 . To apply a quantum algorithm like fac- 
torization to a problem beyond the reach of conventional 
computers, a quantum computer will have to perform 
millions of operations coherently and substantially with- 
out error. Although such control can in principle be 
achieved through quantum error correction |l5| pjj , this 
involves a very high overhead. In contrast, useful quan- 
tum simulations will perhaps require only hundreds of 
operations. 

While it is possible to simulate a quantum system on a 
classical computer, it becomes increasingly difficult as the 
size of the system increases to store the quantum state, 
much less to compute the entire wavefunction evolution. 
For example, a quantum system of 50 spin one-half par- 
ticles occupies a Hilbert space of dimension 2 50 ~ 10 15 . 
This requires ~ 10 15 complex numbers to specify the 
state completely. While intractable using classical com- 
puters, a quantum mechanical device would require only 
50 qubits to store the state of the system. Evolving this 
state vector is also difficult. For the case of local interac- 
tions (all systems obeying special and general relativity), 
Lloyd has suggested an efficient construction of the evo- 
lution operator with small time steps of evolution under 
local interactions [jlCfj . 

This article describes a general scheme for implement- 
ing quantum simulations, and illustrates the flexibility 
of the method with the synthesis of a Hamiltonian not 
normally found in Nature. The challenge that remains 
is to find a sequence of propagators or "gates" that can 
efficiently produce the desired behavior in practical prob- 
lems. 
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II. QUANTUM SIMULATION 

Many of the concepts of quantum simulation are im- 
plicit in the average Hamiltonian theory (AHT) devel- 
oped by Waugh and colleagues to design NMR pulse 
sequences which implement a specific desired effective 
NMR Hamiltonian ||]. While AHT as applied to NMR 
spectroscopy is often directed at obtaining a selectively 
scaled version of the internal Hamiltonian Ti.i„t , the for- 
malism makes clear that other Hamiltonians in the op- 
erator space may be constructed. The significant contri- 
bution to quantum information processing is to articu- 
late the range of propagators U that may be simulated 
given H.i n t and allowed external interactions. Through 
the tenets of quantum information processing it is clear 
that, provided there are coupling pathways (interactions) 
connecting any two identifiable qubits then any Hamilto- 
nian that spans this space may be constructed JT^,^,Q . 
AHT provides a means to quantify the precision of a spe- 
cific simulation and allows a systematic improvement of 
the precision of the implementation. 

A general scheme for quantum simulation utilizing 
some of the results of AHT is based on establishing a 
correspondence between a simulated (model) system, S, 
and a physical (experimental) system, P f2q| . This is 
summarized by the following diagram: 



Ui 



\s(T)} 



\P) 



(2.1) 



\Pi 



The goal is effect the evolution of the simulated system 

\s) — > \ S (T)} using the physical system P, where the 
propagator U = e ~' insT / h 7 and H s is the desired Hamil- 
tonian governing the simulated system. To do this, S is 
related to P by an invertible map 4> which determines 
a correspondence between all the operators and states 
of S and P. In particular, the propagator U maps to 
Vt = </> _1 [/0. The challenge is to implement Vt us- 
ing propagators Vi arising from the available external in- 
teractions with intervening periods of natural evolution 
e -in° p t/h in p SQ that y T = n ie - iW £ ti(T) Fi. If a suffi- 
cient class of simple operations (logic gates) are implc- 
mentable in the physical system, any operator (in partic- 
ular Vt) can be composed of natural evolutions in P and 
external interactions [|10|,E4 26 



. For unitary maps 4>, 
we may write Vt = e~ in P T / h where Ti p = (j^Ti^ can 
be identified with the average Hamiltonian introduced 
by Waugh [HJ29). After \p) \p T ), the final map 

<jr takes \pt) — > \s(T)) thereby effecting the simulation 
s) — » |s(T)). Note that TL S {T) can be a time dependent 
Hamiltonian and that T is treated as a parameter when 



mapped to P. This means that the physical times U{T) 
are parameterized by the simulated time T . 

The desired simulated Hamiltonian may be specified in 
various ways. On the one hand if the simulated Hamilto- 
nian is specified by eigenenergies, translation to a repre- 
sentation in terms of the Pauli matrices o\ may provide 
further physical insight and facilitate implementation via 
geometric algebra techniques p0| using a weakly cou- 
pled spin system. Although no general compiler from a 
Pauli matrix expression is known, the universality of a 
quantum computer implies that an implementation ex- 
ists. Examples of such Hamiltonians include the Balmer 
series of the hydrogen atom and the harmonic oscilla- 
tor. On the other hand, if the simulated Hamiltonian 
is specified in terms of Pauli matrices, translation to an 
eigenstructure representation permits experimental ver- 
ification of its spectral structure. Examples include the 
Ising model TL S — T, i< jaijala J z and the three-body inter- 
action H s = TjirJi23(rl<Tz(jl discussed in this paper. A 
general method for systems with known eigenstructure is 
given in the Appendix. For cases where the eigenstruc- 
ture is not known or the Hamiltonian is not expressed in 
terms of the a\ , simulation is still possible. For example, 
in the treatment of the driven anharmonic oscillator in 
Ref. pa the eigenstructure is not assumed. 



III. IMPLEMENTATION OF EFFECTIVE ctzctzctz 
HAMILTONIAN 



Liquid state NMR quantum computers []3 1| — 133[| are well 
suited for quantum simulations because they have long 
relaxation times (Ti and T2) as well as the flexibility 
of using a variety of molecular samples. In particu- 
lar, the 'scalar' coupling between the nuclear spins, de- 
noted J, may be reduced at will by means of radiofrc- 
quency pulses. Typically spin 1/2 nuclei are used. Thus, 
the kinematics of any 2 N level quantum system could 
be simulated using a given iV-spin molecule. We will 
use a 3-spin system, and illustrate the flexibility of the 
scheme by implementing a non-physical three body in- 
teraction (see also |Q). For a weakly coupled system 
Hint = § PiUW* + ^i<jTr Jij&l&i]- To understand how 
the coupling will behave, we first look at the usual two 
body interaction. 

A scalar two-body coupling propagator of the form 
e -*§ Ji2<? z <? z t ( w fo eie ai are Pauli matrices) will trans- 
form a transverse magnetization, a^, say, into itself and 
an antiphase component: 

(j\ — ^ a\.cosQ + (JyCrl sin#, 

where = irJt. After 8 = ir/2, the antiphase doublet 
<Jy<j1 state is created. An x-phase pulse on spins 1 and 2 
will change this into an antiphase doublet observable on 
spin 2, -a\a-l. 
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Analogously, a three-body coupling propagator of the 
form e~ l ^ Jl23 ' T ' a ''' T ' T (where T is time in the simulated 
system) will transform a transverse magnetization, cr*, 
say, into itself and a component antiphase in the coupled 
spins: 



2 ^123 2 / 

cr — >a^cost 



(3.1) 



After 9 — irJ^T — n/2, the doubly antiphase quartet 
(jy(j1a^ is created. An NMR pulse sequence for imple- 
menting this evolution is straightforward to derive by 
geometric algebra procedures The desired propa- 

gator for the three particle interaction can be expanded 
in terms of the available scalar couplings and free evolu- 
tions: 



3 (-if JllzTa\a\al) 

;(+<f^)g(+ifaJ) e (+<f Jl23T<7>J CT J) fi ( 
;(+«f^) e (+ 4 f CT «) e (+ l f' T ' 



e ^ -4- J ») e (- 

(+if JllzTa\al) x 



;) P (-*f ct x) 



e (+»f ^) e ( +i f CT 55) e (-*§ cr >2) e ^ i 



y> x 



(3.2) 



which implies the following pulse sequence Vt to simulate 
the <j z <j z <j z Hamiltonian: 
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(3.3) 



To generalize, an m-body interaction term can be com- 
posed of a number of two body interaction terms and 
single spin rotations that is linear in m. If not all pairs 
of spins are coupled, as in the case of a linear chain, then 
relay gates must be used |3f| which entails only a poly- 
nomial number of additional operations. 

The three quantum bit NMR system was a room tem- 
perature sample of 13 C labelled Alanine in deuterated 
water. We identify spin 1 as the carbonyl C spin, spin 
2 as the C a spin, and spin 3 as the Cp spin. The scalar 
couplings were J\i = 54.2 Hz, J23 = 35.1 Hz, and J13 = 
1.2 Hz. The 13 C resonance frequency at 9.4 T was 100.6 
MHz, and was detected by an inverse probe. The chemi- 
cal shift difference between spins 1 and 2 was 12,580 Hz, 
and between spins 2 and 3 was 3443 Hz. The proton 
spins were decoupled. Initial states were prepared from 
the thermal equilibrium state with magnetization in all 
three spins by a shaped pulse that excites, for example, 
spins 1 and 3, followed by a magnetic field gradient that 
dephases the magnetization in spins 1 and 3. Then only 
the spin 2 magnetization remains, which may be observed 
by exciting it into a transverse magnetization. Explicitly, 



Here grad refers to a magnetic field gradient, which de- 
stroys the transverse magnetization when viewed as a 
spatial average [ p6[ . 

Figure |l| shows the real and imaginary 13 C spectra 
observed on spin 2 for representative angles 9 or de- 
lay times, confirming that the three particle propagator 
e (-»f Ji23*o- 2 <T z crj transforms the initial state according to 
Equation 3.1. It is clear that the spectra evolve with a 
periodicity T = 2/J123. The simulated time direction 
exhibits only one frequency (positive and negative) since 
there is only one possible non-zero transition energy for 
this system as shown in the Appendix (Eq. AlOj ) . A two- 
dimensional Fourier transform would then directly relate 
the physical eigenenergies to the simulated eigenenergies. 



IV. CONCLUSION 

While multiple quantum coherences have been widely 
used in NMR, as have average Hamiltonian schemes to 
scale down or remove terms from the natural Hamiltonian 
of the physical system ]37|-|39]|, the formalism described 
above adapts these methods to quantum information pro- 
cessing, in particular quantum simulation. Using one 
quantum system to simulate another quantum system 
efficiently is a powerful idea, which may allow the study 
of an otherwise intractable class of problems. Quan- 
tum simulation permits the construction of an arbitrary 
Hamiltonian in the simulated system, as illustrated with 
the construction of the three body interaction a z a z a z 
Hamiltonian. The efficiency, or computational complex- 
ity, of a simulation, however, depends on how difficult 
it is to implement Vt for a given system. Although the 
results given in Q imply that efficient approximate im- 
plementations exist for general cases of physical interest, 
the results discussed in this paper show that more direct, 
exact implementations may also be found in some cases. 
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APPENDIX A: RELATION BETWEEN 
EIGENENERGIES AND PRODUCT OPERATORS 

Representations of Hamiltonians in terms of eigenen- 
ergies may be related to representations in terms of the 
Pauli matrices a\ . Once the arbitrary Hamiltonian is ex- 
pressed in terms of many-body interactions, o],u^. • • ■ cr™, 



,-r 1 _1_ ^-2 , 3 

<7. + (7, + cr. 



0\,, + Cr, 



3 arad 2 
Cr > (7, 



[T/2]J 
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each of these may be broken down in terms of the avail- 
able external and internal Hamiltonians (two body cou- 
plings). Multiple couplings can act at the same time, 
giving a possible increase in efficiency of implementation. 

An arbitrary Hamiltonian, if the eigenstructure is 
known, can be written as 



where E k = E\ E n ~ x ■ ■ ■ E 2 k E\ . Thus we have 



k=0 



The \ipk) are a complete set of orthogonal states and the 
are the energy eigenvalues (not necessarily ordered by 
size). If truncated to the first 2 n levels we have 



2™-l 



n s = th\4>k)(1>k\ 



k=0 



To simulate this with an n-spin system we require the 
map <f> of Eq. [2.1| . One possible map which we adopt is to 
map simulated eigenstates to Zeeman eigenstates: 



\k), 



where \k) is the binary expansion of k. Other mappings 
are possible. The mapping need only connect a basis 
set that spans the 2™ dimensional Hilbert spaces of the 
simulated and physical systems. 

First we resolve the identity in terms of the idempo- 
tents E± = i(l ± a\) where u\ is the Pauli spin matrix 
for spin i: 



1 = {E% + E™){El- 1 
■251 + 



ET 1 ) ■ ■ ■ {E\ + El) 



= E^E^- 1 



E^E^- 1 



... + E n E n-l 
J^El 1 El 1 ' 1 ■ 



£n £n-l 



■E 1 



where the e< take on the values ±1. By identifying 
+ 1^0 and — 1 <-> 1 we may identify the sequence 
En^n-i • • ■ ^2^1 with the binary expansion rjnVn-i ■ ■ ■ mm 
of some integer k, where m take on the values or 1. A 
summary of the relations are: 



e„e„_i • • • e 2 ei <-> mVn-i ■■■mm, 



(Al) 



£i = 1 - 2r]i, 
k = 7 ll 2° + i l2 2 1 

71 

=E^ X - 



Vn2 r 



(A2) 



i=l 

We may then write 



2™-l 



i = J2 E l E Z-\ ■ ■ ■ E i E i E E * 



(A3) 



2"-l 



Us 



n P = E 

k=0 
2™-l 

= E £fc| £ n £ »-l ■ ■ - £ l)( e K £ n~l ■ ■ - e ll 
k=0 
2™-l 

&E e kE e ^ ...E e k. 

k=0 
2™-l 

= E & Ek - ( A4 ) 



fc=0 



In general this will have the form 
H p = (3q 

n n 



j=i j<k=i 

k=0 



(AS) 



where the {a} and {/?} are real numbers. Knowing the 
how do we determine the a k l Substitute a resolution 
of the identity in (|A5|) to get 



2"-l 



2"-l 



k=0 j=0 
j=0 k=0 



(A6) 



Comparing (A6) with (A4) we find that 



2™-l 

This may be written as a matrix equation 
£ = Ma 

where 

M fe , ee (ejW(e^---(e£)< 

Noting that 

2™-l 



(A7) 



MM T = ]T MyAfJ, 

j'=o 

2"-l 

= E M kj M U 



3=0 
2"-l 



fc=0 



= E ( e i) n? ( £ 2 ) vl ■ ■ ■ iA) vl (4)"* • • • 

3=0 
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2"-l 

3=0 
2™-l 

= ^ Ski = 2 n Ski = 2™1, 

3=0 



(A8) 



we see that 



M 



2- n M T . 



Thus inverting the equation (A7) we get the a in terms 
of the energies 



a = 2~ n M T £. 



(A9) 



That is 



3=0 



The matrix M T is the 2™ x 2 n Hadarmard matrix go) 
For example, for two spins, 



M 



1111 

I- 11-1 

II- 1-1 
1-1-1 1 



Therefore, a simulated Hamiltonian specified by £ deter- 
mines the coefficients a in the expansion of the physical 
Hamiltonian with terms that are of the form of many 
body interactions o\a 2 z ■ ■ -a™. For example, the three- 
body interaction Hamiltonian, TL p = ^TrJ^ala^a^. im- 



plies a-j = |7rJi23, and all else is zero. Applying Eq. A9 
to this gives: 



€ = -7rJi23diag(+l, -1, -1, +1, -1, +1, +1, -1) (A10) 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 



[10] 

[11] 
[12] 



P. Benioff, Phys. Rev. Lett. 48, 1581 (1982). 

C. H. Bennett, Int. J. Theor. Phys. 21, 905 (1982). 

D. Deutsch, Proc. R. Soc. Lond. A 400, 97 (1985). 
R. Feynman, Int. J. Theor. Phys. 21, 467 (1982). 
R. Landauer, Phys. Lett. A 217, 188 (1996). 

S. Lloyd, Science 261, 1569 (1993). 

S. Lloyd, Phys. Rev. Lett. 75, 4714 (1995). 

P. W. Shor, SI AM J. Comput. 26, 1484 (1997). 

J. F. Traub, G. W. Wasilkowski, H. Wozniakowski, 

Information-Based Complexity Academic Press, New 

York (1988). 

S. Lloyd, Science 273, 1073 (1996). 

S. Weisner, xxx. lanl.gov/quant-ph/9603028 (1996). 

C. Zalka, Proc. Roy. Soc. Lndn A 454, 313 (1998). 



[13] D. S. Abrams and S. Lloyd, Phys. Rev. Lett. 79 2586 

(1997) . 

[14] D. A. Lidar and O. Biham, Phys. Rev. E 56, 3661 (1997). 
[15] E. Knill, R. Laflamme and W. Zurek, Science , 279, 342 

(1998) . 

[16] E. Knill and R. Laflamme, Phys. Rev. A 55, 900 (1997). 
[17] D. Cory, M. Price, W. Maas, E. Knill, R. Laflamme, W. 

Zurek, T. Havel, and S. Somaroo, Phys. Rev. Lett. 81, 

2152 (1998). 
[18] A. M. Steane, Nature 399, 124 (1999). 
[19] S. L. Braunstein, Nature 394, 47 (1998). 
[20] A. R. Calderbank , E. M. Rains, P. W. Shor, N. J. A. 

Sloane, Phys. Rev. Lett. 78 405 (1997). 
[21] R. Laflamme, C. Miquel, J. Pablo Paz, W. H. Zurek, 

Phys. Rev. Lett. 77 198 (1996). 
[22] U. Haeberlen and J. Waugh, Phys. Rev. 175, 453 (1968). 
[23] Barenco A, Bennett CH, Cleve R, DiVincenzo DP, Mar- 

golus N, Shor P, Sleator T, Smolin JA, Weinfurter H, 

Phys. Rev. A 52, 3457 (1995). 
[24] D.P. DiVincenzo, Science 270, 255 (1995). 
[25] S. Somaroo, CH. Tseng, T. F. Havel, R. Laflamme, D.G. 

Cory, Phys. Rev. Lett. 82, 5381 (1999). 
[26] A. Barenco, C. Bennett, R. Cleve, D. DiVincenzo, N. 

Margolus, P. Shor, T. Sleator, J. Smolin and H. Wein- 
furter, Phys. Rev. A. 52, 3457 (1995). 
[27] D. P. DiVincenzo, Phys. Rev. A 51, 1015 (1995). 
[28] D. P. Divincenzo and P. W. Shor, Phys. Rev. Lett. 77, 

3260 (1996). 

[29] J. Ellett and J. S. Waugh, J. Chem. Phys. 51, 2581 
(1969). 

[30] S. Somaroo, D. Cory and T. Havel, Phys. Lett. A 240, 
1 (1998). 

[31] D. Cory, A. Fahmy and T. Havel, Proc. Natl. Acad. Sci. 

94, 1634-1639 (1997). 
[32] N. Gershenfeld and I. Chuang, Science 275, 350-356 

(1997) . 

[33] A. Steane, Phys. Rev. Lett. 77, 793 (1996). 

[34] R. Bruschweiler, J. Chem. Phys. 105 , 6164 (1996). 

[35] J. Evans, Biomolecular NMR Spectroscopy, p. 64, Oxford 

University Press, Oxford, (1995). 
[36] D. G. Cory, M. D. Price, T. F. Havel, Physica D 120, 

82 (1998). 

[37] W. S. Warren, J. Chem. Phys. 93, 858 (1990). 

[38] M. Eden and M. H. Levitt, Chem. Phys. Lett. 293 173 

(1998) . 

[39] M. H. Levitt and R. R. Ernst, Journal of Chem. Phys. 

83 3297(1985). 

[40] J. Hadamard, Bull. Sciences Math. 17, 240 (1893). 



5 



Real 



Imaginary 





FIG. 1. NMR spectra from C a spin of 1J C-labelled Alanine 
demonstrating a quantum simulation of the o z o z o z Hamilto- 
nian as implemented by the pulse sequence (3.3). The C a 
resonance is split twice by the couplings to the other two car- 
bon nuclei, resulting in four lines. As a function of the angle 
(or evol ution time) the spectra exhibit the periodicity given 
by (3.1): 6 = n/2, doubly antiphase; = tt, in phase nega- 
tive; 6 — 3n/2, doubly antiphase; 9 — 2n, back to in phase 
positive. 
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